Search for a diffuse flux of high-energy with the ANTARES 

neutrino telescope 



J.A. Aguilar 3 , I. Al Samarai b , A. Albert , M. Andre d , M. Anghinolfi 6 , G. Anton f , S. Anvar g , M. 
Ardid h , A.C. Assis Jesus', T. Astraatmadja u , J-J. Aubert b , R. Auer f , B. Baret j , S. Basa k , M. 
Bazzotti 1 ™, V. Bertin b , S. Biagi 1 " 1 *, C. Bigongiari 3 , C. Bogazzi', M. Bou-Cabo h , M.C. 
Bouwhuis 1 , A.M. Brown b , J. Brunner b ' 2 , J. Busto b , F. Camarena h , A. Capone" , C.Carloganu p , 

G. Carminati l m , J. Carr b , S. Cecchini l q , Ph. Charvis 1 , T. Chiarusi 1 , M. Circella 5 , R. Coniglione v , 

H. Costantini 6 , N. Cottini', R Coyle b , C. Curtil b , M.R Decowski', I. Dekeyser 11 , A. Deschamps 1 , 
C. DonzaudJ ™, D. Dornic 3 , Q. Dorosti x , D. Drouhin , T. Eberl f , U. Emanuele 3 , J-R Ernenwein b , 
S. Escoffier b , F. Fehr f , V. Flaminio y z , F. Folger f , U. Fritsch f , J-L. Fuda", S. Galata b , P. Gay p , G. 
Giacomelli lm , J.R Gomez-Gonzalez 3 , K. Graff, G. Guillard 33 , G. Halladjian b , G. Hallewell b , H. 

van Haren 3b , A.J. Heijboer', Y. Hello 1 , J.J. Hernandez-Rey 3 , B. Herold f , J. Ho61 f , C.C. Hsu', 
M. de Jong' 1 , M. Kadler 3C , N. Kalantar-Nayestanaki x , O. Kalekin f , A. Kappes f , U. Katz f , P. 
Kooijman 1,ad ' 3e , C. Kopper', A. KouchnerJ, V. Kulikovskiy af e , R. Lahmann', P. Lamare 8 , G. 
Larosa h , D. Lefevre", G. Lim Ue , D. Lo Presti 3g , H. Loehner x , S. Loucatos 1 , F. Lucarelli"' , S. 

Mangano 3 , M. Marcelin\ A. Margiotta 1 '" 1 , J.A. Martinez-Mora h , A. Mazure k , A. Meli f , T. 
Montaruli 8 ' 311 , M. Morganti y ' z , L. Moscoso t,j , H. Motz f , C. Naumann 1 , M. Neff*, D. Palioselitis', 
G.E.Pavalas 3 ', P. Payre b , J. Petrovic', P. Piattelli v , N. Picot-Clemente b , C. Picq 1 , V. Popa 3 ', T. 

Pradier 33 , E. Presani', C. Racca c , C. Reed', G. Riccobene v , C. Richai-dt f , K. Roensch f , A. 
Rostovtsev 3j , M. Rujoiu 3 ', G.V. Russo 3g , F. Salesa 3 , P. Sapienza v , F. Sch6ck f , J-P. Schuller 1 , R. 

Shanidze f , F. Simeone"' , A. Spies*, M. Spurio 1 ™'*, J.J.M. Steijger', Th. Stolarczyk 1 , M. 
Taiuti e ' 3k , C. Tamburini 11 , L. Tasca k , S. Toscano 3 , B. Vallage 1 , V. Van Elewyck, G. Vannoni 1 , M. 
Vecchi n b , P. Vernin 1 , G. Wijnker 1 , E. de Wolf- 36 , H. Yepes 3 , D. Zaborov aj , J.D. Zornoza 3 , 

J. Zufiiga 3 

a IFIC - Instituto de Fisica Corpuscular, Edificios Investigation de Paterna, CSIC - Universitat de Valencia, Apdo. de Correos 22085, 

46071 Valencia, Spain 
b CPPM, Aix-Marseille Universite, CNRS/IN2P3, Marseille, France 
L 'GRPHE - Institut universitaire de technologic de Colmar, 34 rue du Grillenbreit BP 50568 - 68008 Colmar, France 
Technical University of Catalonia, Laboratory of Applied Bioacouslics.Rambla Exposicio, 08800 Vilanova i la Geltru,Barcelona, Spain 
e INFN - Sezione di Genova, Via Dodecaneso 33, 16146 Genova, Italy 
^Friedrich-Alexander-Universitdt Erlangen-Niirnberg, Erlangen Centre for Astroparticle Physics, Erwin-Rommel~Str. 1, 91058 

Erlangen, Germany 

^Direction des Sciences de la Matiere - Institut de recherche sur les lois fondamentales de VUnivers - Service d'Electronique des 
Detecteurs et d'Informatique, CEA Saclay, 91191 Gif-sur-Yvette Cedex, France 
h Institut d'lnvestigacio per a la Gestio Integrada de Zones Costaneres (IGIC) * Universitat Politecnica de Valencia. C/ Paranimf 1. , 

46730 Gandia, Spain. 
'Nikhef Science Park, Amsterdam, The Netherlands 
1 APC - Laboratoire AstroParticule et Cosmologie, UMR 7164 ( CNRS, Universite Paris 7 Diderot, CEA, Observatoire de Paris) 10, rue 
Alice Domon et Leonie Duquet 75205 Paris Cedex 13, France 
k LAM ~ Laboratoire d'Aslrophysique de Marseille, Pole de VEtoile Site de Chateau-Gombert, rue Frederic Joliot-Curie 38, 13388 

Marseille Cedex 13, France 
'INFN - Sezione di Bologna, Viale Berti Pichat 6/2, 40127 Bologna, Italy 
'"Dipartimento di Fisica dell'Universita, Viale Berti Pichat 6/2, 40127 Bologna, Italy 
" Dipartimento di Fisica dell'Universita La Sapienza, Pie Aldo Mom 2, 00185 Roma, Italy 
"INFN -Sezione di Roma, P. le Aldo Mow 2, 00185 Roma, Italy 
p Laboratoire de Physique Corpusculaire, 1N2P3-CNRS, Universite Blaise Pascal, Clermont-Ferrand, France 
1INAF-IASF via P. Gobetti 101, 40129 Bologna, Italy 
r Geoazur - Universite de Nice Sophia-Antipolis, CNRS/1NSU, IRD, Observatoire de la Cote d'Azur and Universite Pierre et Marie Curie, 

BP 48, 06235 Villefranche-sur-mer, France 
S INFN - Sezione di Bari, Via E. Orabona 4, 70126 Bari, Italy 
'Direction des Sciences de la Matiere ■ Institut de recherche sur les lois fondamentales de VUnivers - Service de Physique des Particules, 

CEA Saclay, 91191 Gif-sur-Yvette Cedex, France 
"COM - Centre d'Oceanologie de Marseille, CNRS/INSU et Universite de la Mediterranee, 163 Avenue de Luminy, Case 901, 13288 

Marseille Cedex 9, France 



1 



'INFN - Labomtori Nazionali del Sud (LNS), Via S. Sofia 62, 95123 Catania, Italy 
w Univ Paris-Sud, 91405 Orsay Cedex, France 
x Kernfysisch Versneller Instituut (KV1), University ofGroningen, Zernikelaan 25, 9747 AA Groningen, The Netherlands 
y Dipartimento di Fisica delVUniversita, Largo B. Pontecorvo 3, 56127 Pisa, Italy 
Z INFN - Sezione di Pisa, Largo B. Pontecorvo 3, 56127 Pisa, Italy 
a IPHC-Institut Pluridisciplinaire Hubert Curien - Universile de Strasbourg et CNRS/IN2P3 23 rue du Loess, BP 28, 67037 Strasbourg 

Cedex 2, France 

"''Royal Netherlands Institute for Sea Research (NIOZ), Landsdiep 4,1797 SZ 't Horntje (Texel), The Netherlands 
ac Dr Remeis Sternwarte Bamberg, Slernwartstrasse 7, Bamberg, Germany 
a<1 Universiteit Utrecht, Faculteit Betawetenschappen, Princetonplein 5, 3584 CC Utrecht, The Netherlands 
ae Universiteit van Amsterdam, Instituut voor Hoge-Energie Fysika, Science Park 105, 1098 XG Amsterdam, The Netherlands 
a f Moscow State University, Skobellsyn Institute of Nuclear Physics, Leninskie gory, 119991 Moscow, Russia 
ag Dipartimento di Fisica ed Astronomia delVUniversita, Viale Andrea Doria 6, 95125 Catania, Italy 
ah University of Wisconsin - Madison, 53715, WI, USA 
al Institute for Space Sciences, R-77125 Bucharest, Magurele, Romania 
aJ ITEP - Institute for Theoretical and Experimental Physics, B. Cheremushkinskaya 25, 117218 Moscow, Russia 
ak Dipartimento di Fisica delVUniversita, Via Dodecaneso 33, 16146 Genova, Italy 



Abstract 

A search for a diffuse flux of astrophysical muon neutrinos, using data collected by the ANTARES 
neutrino telescope is presented. A (0.83 x 2n) sr sky was monitored for a total of 334 days of 
equivalent live time. The searched signal corresponds to an excess of events, produced by astro- 
physical sources, over the expected atmospheric neutrino background. The observed number of 
events is found compatible with the background expectation. Assuming an E~ 2 flux spectrum, 
a 90% c.l. upper limit on the diffuse flux of E 2 <&c, m = 5.3 x 10~ 8 GeV cnT 2 s _1 sr~' in the 
energy range 20 TeV - 2.5 PeV is obtained. Other signal models with different energy spectra 
are also tested and some rejected. 
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1. Introduction 



This letter presents a search for a diffuse flux of high energy muon neutrinos from astrophys- 
ical sources with the ANTARES neutrino telescope. The construction of the deep sea ANTARES 
detector was completed in May 2008 with the connection of its twelfth detector line. The tele- 
scope is located 42 km off the southern coast of France, near Toulon, at a maximum depth of 
2475 m. 

The prediction of the diffuse neutrino flux from unresolved astrophysical sources is based on 
cosmic ray (CR) and y-ray observations. Both electrons (leptonic models) HHH and protons or 
nuclei (hadronic models) |3[] can be accelerated in astrophysical processes. In the framework of 
hadronic models the energy escaping from the sources is distributed between CRs, y-rays and 
neutrinos. Upper bounds for the neutrino diffuse flux are derived from the observation of the 
diffuse fluxes of y-rays and ultra high energy CRs taking into account the production kinematics, 
the opacity of the source to neutrons and the effect of propagation in the Universe. There are two 
relevant predictions: 

- The Waxman-Bahcall (W&B) upper bound [4] uses the CR observations at Eqr ~ 10 19 eV 
(Ec R Q>cr ~ 10~ 8 GeV cm _2 s _1 sr _1 ) to constrain the diffuse flux per neutrino flavour (here and in 
the following the symbol v represents the sum of plus v^,): 

£ 2 O v < 4.5/2 x 10~ 8 GeV cnr 2 sr _1 s _1 (1) 

(the factor 1/2 is added to take into account neutrino oscillations). This value represents a bench- 
mark flux for neutrino telescopes. 

- The Mannheim-Protheroe-Rachen (MPR) upper bound $^ is derived using as constraints the 
observed CR fluxes over the range from 10 5 to 10 9 GeV and y-ray diffuse fluxes. In the case 
of sources opaque to neutrons, the limit is E 2 O v < 2 x 10~ 6 ( GeV cm _2 sr _1 s _1 ); in the case 
of sources transparent to neutrons, the limit decreases from the value for opaque sources at 
E v ~ 10 6 GeV to the value of Eq. [Tjat E v ~ 10 9 GeV. 

The detection of high energy cosmic neutrinos is not background free. Showers induced by 
interactions of CRs with the Earth's atmosphere give rise to atmospheric muons and atmospheric 
neutrinos. Atmospheric neutrinos that have traversed the Earth and have been detected in the 
neutrino telescope, are an irreducible background for the study of cosmic neutrinos. As the 
spectrum of cosmic neutrinos is expected to be harder (oc E~ 2 ) than that of atmospheric neutrinos, 
a way to distinguish the cosmic diffuse flux is to search for an excess of high energy events in 
the measured energy spectrum. 

The relevant characteristics of the ANTARES detector are presented in Sec. 2. The rejection 
of the atmospheric muon background and an estimator of the muon energy are discussed in Sec. 
3. This estimator is used to discriminate high energy neutrino candidates from the bulk of lower 
energy atmospheric neutrinos. The results are presented and discussed in Sec. 4 and Sec. 5. 

2. reconstruction in the ANTARES detector 

The ANTARES detector is a three-dimensional array of photomultiplier tubes (PMTs) dis- 
tributed along twelve lines [6]. Each line comprises 25 storeys, spaced vertically by 14.5 m, 
with each storey containing three optical modules (OMs) 071] and a local control module for the 
corresponding electronics [8]. The OMs (885 in total) are arranged with the axes of the PMTs 
oriented 45° below the horizontal. The lines are anchored on the seabed at distances of about 70 
m from each other and tensioned by a buoy at the top of each line. 
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Muon neutrinos are detected via charged current interactions: +N — > fX+X. The challenge 
of measuring muon neutrinos consists of reconstructing the trajectory using the arrival times 
and the amplitudes of the Cherenkov light signal detected by the OMs, and of estimating the 
energy. The track reconstruction algorithm ]9J is based on a likelihood fit that uses a detailed 
parametrization of the probability density function for the photon arrival times taking into ac- 
count the delayed photons. The output is: the track position and direction; the information on the 
number of hits (Ni,i t ) used for the reconstruction; a quality parameter A. A is determined from 
the likelihood and the number of compatible solutions found by the algorithm and can be used 
to reject badly reconstructed events. Without any cut on A, the fraction of atmospheric muon 
events that are reconstructed as upward-going is ~ 2% (see Table [TJ. The appropriate value of 
the A variable cut for this analysis is discussed in Sec. 3. Monte Carlo (MC) simulations show 
that the ANTARES detector achieves a median angular resolution for muon neutrinos better than 
0.3° for £ v > lOTeV. 

Muon energy losses are due to several processes ifioll and can be parametrized as: 

dE M Idx = aiE^) + f3(E M ) ■ E M , (2) 

where aiE^) is an almost constant term that accounts for ionisation, and PiE^) takes into account 
the radiative losses that dominate for Ep >0.5 TeV. Particles above the Cherenkov threshold 
produce a coherent radiation emitted in a Cherenkov cone with a characteristic angle 6>c - 43° in 
water. Photons emitted at the Cherenkov angle, arriving at the OMs without being scattered, are 
referred to as direct photons. The differences between the calculated and the measured arrival 
time (time residuals) of direct photons follow a nearly Gaussian distribution of few ns width, due 
to the chromatic dispersion in the sea water and to the transit time spread of the PMTs. 

For high muon energies (E M > 1 TeV), the contribution of the energy losses due to radiative 
processes increases linearly with the muon energy and the resulting electromagnetic showers 
produce additional light. 

Scattered Cherenkov radiation or photons originating from secondary electromagnetic show- 
ers arriving on the OMs (denoted from now on as delayed photons), are delayed with respect to 



the direct photons, with arrival time differences up to hundreds of ns 111 111 . As a consequence, the 



percentage of delayed photons with respect to direct photons increases with the muon energy. 



The PMT signal is processed by two ASIC chips (the Analogue Ring Sampler, ARS [12]) 
which digitize the time and the amplitude of the signal (the hit). They are operated in a token 
ring scheme. If the signal crosses a preset threshold, typically 0.3 photo-electrons, the first ARS 
integrates the pulse within a window of 25 ns and then hands over to the second chip with a 
dead time of 15 ns. If triggered, the second chip provides a second hit with a further integration 
window of 25 ns. After digitization, each chip has a dead time of typically 250 ns. After this 
dead time, a third and fourth hit can also be present. 



2.1. The Monte Carlo simulations 

The simulation chain lf]~3L [l4ll comprises the generation of Cherenkov light, the inclusion of 
the optical background caused by bioluminescence and radioactive isotopes present in sea water, 
and the digitization of the PMT signals. Upgoing muon neutrinos and downgoing atmospheric 
muons have been simulated and stored in the same format used for data. 

Signal and atmospheric neutrinos. MC muon neutrino events have been generated in the en- 
ergy range 10 < E v < 10 8 GeV and zenith angle between 0° < 9 < 90° (upgoing events). 
The same MC sample can be differently weighted to reproduce the "conventional" atmospheric 
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neutrinos from charged meson decay (Bartol) [15] (<1> V oc E V 3J at high energies), the "prompt" 
neutrinos and the theoretical astrophysical signal (O v oc E~ 2 ). A test spectrum: 

El <£ v = 1.0 x KT 7 GeVcnT 2 s" 1 sr _1 . (3) 

is used to simulate the diffuse flux signal. The normalization of this test flux is irrelevant when 
defining cuts, optimizing procedures, and calculating the sensitivity. 

Above 10 TeV, the semi-leptonic decay of short-lived charmed particles D — > K + fi + v M 
becomes a significant source of atmospheric "prompt leptons". The lack of precise information 
on high-energy charm production in hadron-nucleus collisions leads to a great uncertainty (up 
to four orders of magnitude) in the estimate of the leptonic flux above 100 TeV. The models 
considered in I1 1611 were used, in particular the Recombination Quark Parton Model (RQPM) 
which gives the largest prompt contribution. 

Atmospheric muons. Atmospheric muons reconstructed as upgoing are the main background 
for a neutrino signal and their rejection is a crucial point in this analysis. Atmospheric muon 
samples have been simulated with the MUPAGE package II 1711 . In addition to one month of 



equivalent live time with a total energy Et > 1 GeV 111811 . a dedicated one year of equivalent live 
time with Ej > 1 TeV and multiplicity m — 1 H- 1000 was generated. The total energy Ej is the 
sum of the energy of the individual muons in the bundle. Triggered ANTARES events mainly 
consist of multiple muons originating in the same primary CR interaction. For the ANTARES 
detector the background contribution from muon events originating from independent showers is 
negligible. 

Simulation of the detector. In the simulation of the digitized signal, the main features of the 
PMTs and of the ARSs are taken into account. The simulated photons arriving on each PMT 
are used to determine the charge of the analogue pulse; the charge of consecutive pulses are 
added during the 25 ns integration time. The hit time is determined by the arrival time of the 
first photon. To simulate the noise in the apparatus, background hits generated according to the 
distribution of a typical data run are added. The status of each of the 885 OMs in this particular 
run is also reproduced. The OM simulation also includes the probability of a detected hit giving 
rise to an afterpulse in the PMT. This probability was measured in the laboratory 111911 and was 
confirmed with deep-sea data. 



3. Event selection and background rejection 

The data were collected during the period from December 2007 to December 2009 with 9, 
10 and 12 active line configurations. The runs were selected according to a set of data-quality 
criteria described in Hill : in particular a baseline rate < 120 kHz and a burst fraction < 40%. 
A total of 3076 runs satisfy the conditions. The total live time is 334 days: 70 days with 12 
lines, 128 days with 10 lines and 136 days with 9 lines. In the detector simulation three different 
configurations are taken into account, based on the number of active lines. For each detector 
geometry, a typical run is selected to reproduce on average the conditions and the background in 
the data. 

3.1. Rejection of atmospheric muons 

The ANTARES trigger rate, which is dominated by atmospheric muons, is a few Hz. The 
reconstruction algorithm ^ results in approximately 5% of triggered downgoing muons to be 
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HAtm 


VAtm 




Data 


Reco 


2.2 x 10 8 


7.11 X 10 3 


106 


2.5 x 10 s 


Upgoing 


4.8 x 10 6 


5.50 x 10 3 


80 


5.2 x 10 6 


1 "-level 


9.1 x 10 3 


142 


24 


1.0 x 10 4 


2'"'-level 





116 


20 


134 



Table 1: Expected event number in 334 days of equivalent live time for the three MC samples (atmospheric 
muons, atmospheric neutrinos (Bartol+RQPM), astrophysical signal from Eq. [3] and data. Reco: at the 
reconstruction level; Upgoing: reconstructed as upgoing; l"-level: after the first-level cuts; 2'"'-level: after 
the second-level cut. The number of events in the data set is known only at the end, after the un-blinding 
procedure, see Sec. 4. 



mis-reconstructed as upgoing. This contamination can be readily reduced by applying require- 
ments on the geometry of the event and on the track reconstruction quality parameter A. For 
simulated upgoing atmospheric neutrino events the A distribution has a maximum around -4.5 
and 95% of the events have A > -5.5. Two steps are used to remove the contamination of 
mis-reconstructed atmospheric muons from the final sample. 

First-level cuts. Selection of (z) upgoing particles with reconstructed zenith angle 9 rec < 80° 
(corresponding to 0.83 x27r sr); (z'z) A > -6; (z'zz) Nut > 60; (z'v) reconstruction with at least two 
lines. The first-level cuts remove all MC atmospheric muons with Ej < 1 TeV and reduce the 
rate of mis-reconstructed events by almost 3 orders of magnitude, as indicated in Table Q] 
Second-level cut. The remaining mis-reconstructed atmospheric muons have a quality parameter 
A which on average decreases with increasing Nhn- Values of a cut parameter (A*) are obtained 
in 10 different intervals of Nhu, in order to reduce the expected rate of mis-reconstructed events 
to less than 0.1 -r 0.3 events/year in each interval. A parametrization of the values of A* as a 
function of Nhu is: 

f -4.59 - 5.88 ■ W- 3 N hit for N m < 172 

\ -5.60 for N hit > 172 ( ' 

Removing all events with A < A*, the atmospheric muons are completely suppressed (last row 
of Table [T). Independent MC atmospheric muon simulations using CORSIKA (see details in 



1 1411) confirm that the maximum contamination in the final sample is less than 1 event/year. The 
effects of the first- and second-level cuts on signal and atmospheric neutrinos are also given in 
Tabled] 



3.2. The energy estimator 

To separate atmospheric and astrophysical neutrinos, an original energy estimator is defined, 

which is based on hit repetitions in the OMs due to the different arrival time of direct and delayed 

photons. The number of repetitions R, for the z'-th OM is defined as the number of hits in the same 

OM within 500 ns from the earliest hit selected by the reconstruction algorithm. In most cases, 

y r 

Rj =1 or 2. The mean number of repetitions in the event is defined as R — where Nom is the 
number of OMs in which hits used by the tracking algorithm are present. After the second-level 
cut, R is linearly correlated with the log of the true muon energy E true in the range from 10 TeV 
to 1 PeV, see Fig. Q] R slightly saturates after 1 PeV. The distribution of log(E rec /E, me ) has a 
HWHM=0.4 when R is used as an estimator of the muon energy E rec . This energy estimator 
is robust because it does not depend on the number of active OMs and on non-linear effects on 
charge integration. 
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2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7 

log io (E^/GeV) 

Figure 1: Mean number of repetitions as a function of the true neutrino-induced muon energy for events 
passing the second-level cut. 

Atmospheric muons are used to check the agreement between data and MC for the R variable. 
The first-level cuts (Sec. 13. U are applied both to data and MC, except that tracks reconstructed as 
downgoing (8 rec > 90°) are selected. In the data set, 1 .37 x 10 7 events are present, and 1 .22 x 10 7 
in the simulation for the corresponding live time. The distribution of the R variable is shown in 
Fig. |2]a) for a subset of 20 days live time of the 12 line data. A second comparison in Fig. |2]b) 
uses those atmospheric muons that survived the first-level cuts in the same data set and are mis- 
reconstructed as upgoing (the true upgoing atmospheric neutrinos are about 1.5% of the total, see 
Table[T]i. The MC curve a) was normalized to the data by a factor 1.12, and b) by a factor 1.15. 
These factors are well within the overall systematic uncertainties on the atmospheric muon flux, 
and the relative difference is accounted for by the uncertainty on the OMs angular acceptance 
Q. 

3.3. Signal/atmospheric v M background discrimination 

The separation of the diffuse flux signal from the atmospheric background is performed 
by a cut on the R variable. In order to avoid any bias, a blinding procedure on MC events is 
applied, without using information from the data. The numbers of expected events for signal (n s ) 
and background («/,) are computed as a function of R to find the optimal cut value of R. Later 
the number of observed data events (n i, s ) are revealed (un-blinding procedure) and compared 
with the expected background for the selected region of R. If this number is compatible with the 
background, the upper limit for the flux at a 90% confidence level (c.l.) is calculated using the 
Feldman-Cousins method 1 20] . 

Simulated atmospheric neutrino events are used also to calculate the "average upper limit" 
that would be observed by an ensemble of hypothetical experiments with no true signal (n s = 0) 
and expected background «/,. Taking into account all the possible fluctuations for the estimated 
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background, weighted according to their Poisson probability of occurrence, the average upper 
limit is: 

Jho%(- n b)= Y\ fi90%(n obs ,n b ) ^ e~ n ». (5) 

n ( n obs)- 

n„b s =v 



The best average upper limit is obtained with the cut on the energy estimator that minimizes the 
so-called M( 
upper limit: 



so-called Model Rejection Factor [21], MRF = an( i hence minimizes the average flux 



O 90% =O v -^^«^=avMRF. (6) 

n s 

The value of R which minimizes the MRF function in Eq. [6] is used as the discriminator 
between low energy events, dominated by the atmospheric neutrinos, and high energy events, 
where the signal could exceed the background. 

The method relies on knowledge of the number of background events expected for a given 
period of data. The cumulative distributions of the R variable are computed for atmospheric neu- 
trino background and diffuse flux signal for the three discussed configurations of the ANTARES 
detector and the corresponding live times. For the atmospheric neutrino background, the con- 
ventional flux and the prompt models are considered separately. Fig. [3] shows the cumulative 
distributions of the R variable for signal and background neutrinos (Bartol+RQPM). Using these 
cumulative distributions, the MRF is calculated as a function of R; the minimum (MRF=0.65) 
is found for R = 1.31. Assuming the Bartol (Bartol+RQPM) atmospheric fluxes, 8.7 (10.7) 
background events and 10.8 signal events (assuming the test flux of Eq. [3) are expected for 
R > 1.31. Fig. |4] shows the energy spectra for signal and background neutrino events before 




Figure 3: Cumulative distributions of the R variable for simulated diffuse flux signal (Eq. |3) and atmospheric 
neutrino events (including the prompt from the RQPM model). 



and after the cut R > 1.31. The central 90% of the signal is found in the neutrino energy range 
20 TeV < E v < 2.5 PeV. 



Flux model 


N <U i 


D/N <L31 


N>,.3, 


Kim 


Bartol 


104.0 


1.20 


8.7 


10.4 


Bartol+ RQPM 


105.2 


1.19 


10.7 


12.7 



Table 2: Number of expected atmospheric neutrino events in the two intervals of R. Atmospheric v,, accord- 
ing to the Bartol flux and the Bartol flux plus the prompt contribution from the RQPM model are considered 
separately. N<[ 31 is for R < 1.31, where D=125 data events are observed. N sl 31 (N* , 31 ) is the number of 
expected background events for R > 1.31, without (with) the normalization by the factor D/N<i3i. 



4. Data un-blinding and results 

Events surviving the second-level cut are upgoing neutrino candidates. Fig. [5] shows the 
distribution of the neutrino candidates as a function of R, compared with that given by the atmo- 
spheric neutrino MC. At this stage, only the 125 events with R < 1.31 are un-blinded. The events 
with R > 1.31 in Fig. [5]are revealed only after the un-blinding of the data samples. The number 
of expected events is lower by ~ 20% with respect to the detected events (D). This discrepancy is 
well within the systematic uncertainties of the absolute neutrino flux at these energies (25-30%) 

Table |2] shows the number of expected MC events N>i .31 and N* j 31 = N>i .31 • D/N<i 3i 
both for the conventional Bartol and Bartol+RQPM fluxes. Most prompt models give negligible 
contribution (the average over all considered models gives 0.3 events), the RQPM model predicts 
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Atmsv, no R cut 




2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7 

log io (E v /GeV) 



Figure 4: Signal and background neutrino energy spectra as a function of the true neutrino energy after the 
second-level cut with and without the requirement R > 1.31. The energy region containing 90% of the 
astrophysical neutrino signal is indicated. 



the largest contribution of 2.0 additional events with respect to the conventional Bartol flux. After 
data/MC normalization in the R < 1.31 region, the number of expected background events for 
R > 1.31 from a combined model of Bartol flux plus the average contribution from prompt 
models is 10.7 events. 

A reasonable agreement between data and MC for the R distribution both for atmospheric 
muons (c.f. Fig. |2]i and for atmospheric neutrinos in the test region R <1.31 (c.f. Fig. |5]l is 
found. Consequently the data was un-blinded for the signal region R > 1.31 and 9 high-energy 
neutrino candidates are found. 

Systematic uncertainties on the expected number of background events in the high energy 
region (R > 1.31) include: (z) the contribution of prompt neutrinos, estimated as *H events. In 
the following, the largest value is conservatively used, (z'z") The uncertainties from the neutrino 
flux from charged meson decay as a function of the energy. By changing the atmospheric neutrino 
spectral index by +0.1, both below and above ~10 TeV (when the conventional neutrino flux 
has spectral index one power steeper than that of the primary CR below and after the knee, 
respectively), the relative number of events for R > 1.31 changes at most by +1.1, keeping in 
the region R < 1.31 the number of MC events equal to the number of data. The migration from 
the Bartol to the Honda MC [22] produces a smaller effect. The uncertainties on the detector 
efficiency (including the angular acceptance of the optical module 01411 . water absorption and 
scattering length, trigger simulation and the effect of PMT afterpulses) amount to 5% after the 
normalization to the observed atmospheric background in the test region. 

The number of observed events is compatible with the number of expected background 
events. The 90% c.l. upper limit on the number of signal events fXgo%(nb) for "fc = 10.7 ± 2 back- 
ground events and n b s = 9 observed events including the systematic uncertainties is computed 
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CM 




Figure 5: Distribution of the R parameter for the 134 neutrino candidates in the 334 days of equivalent live 
time. Points represent data, the filled histogram is the atmospheric neutrino MC (Bartol model only). The 
dashed line represents the maximum contribution (RQPM) of "prompt" neutrinos. The MC predictions are 
not normalized to the data. The signal at the level of the upper limit (Eq. [7j( is shown as a full line. The cut 
at R = 1.31 is indicated as a vertical line. 



with the method of B23I1 . The value /U9o%(«fe) = 5.7 is obtained. The profile likelihood method 



H2411 gives similar results. The corresponding flux upper limit is given by <&9o% = O v ■ fi9o%/n s : 



E 2 ® m = 5.3 x 1(T 8 GeV cm" 2 s" 1 sr" 1 (7) 

(our expected sensitivity is 7.0 x 10~ 8 GeV crrT 2 s _1 sr~'). This limit holds for the energy range 
between 20 TeV to 2.5 PeV, as shown in Fig. [4] The result is compared with other measured flux 
upper limits in Fig. 

A number of models predict cosmic neutrino fluxes with a spectral shape different from E~ 2 . 
For each model a cut value R* is optimized following the procedure in Sec. 3.3. Table[3]gives the 
results for the models tested; the value of R* ; the number N mo d of signal events for R > R*; the 
energy interval where 90% of the signal is expected; the ratio between figo% (computed according 
to 112011 ) and N moi j. A value of ^o%/N mo d < 1 indicates that the theoretical model is inconsistent 
with the experimental result at the 90% c.l. In all cases (except for B33IO . our results improve 



upon those obtained in Il27l.l28l. l2' 



3 Charged current v r interaction can contribute via r~ — » i-Cvjv^ (and similarly the v T ) by less than ~ 10% both for 
signal and background. For the background, the v T contribution is almost completely absorbed by the uncertainty on the 
overall normalization, while it is neglected in the signal. 
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Figure 6: The ANTARES 90% c.l. upper limit for a E~ 2 diffuse high energy v,, + flux obtained in this work, compared 
with the limits from other experiments. The Frejus [25], MACRO [26], Amanda-II 2000-03 [27] limits refer to v„ + v„. 
The Baikal [ 28] and Amanda-II UHE 2000-02 [29] refer to neutrinos and antineutrinos of all-flavours, and are divided by 
3. For reference, the W&B [ 4] and the MPR| 5] upper bounds for transparent sources are also shown. They are divided 
by two, to take into account neutrino oscillations. The grey band represents the expected variation of the atmospheric 
flux: the minimum is the Bartol flux from the vertical direction; the maximum the Bartol+RQPM flux from the horizontal 
direction. The central line is averaged over all directions. 



5. Conclusions 

A search for a diffuse flux of high energy muon neutrinos from astrophysical sources with 
the data from 334 days of live time of the ANTARES neutrino telescope is presented. A robust 
energy estimator, based on the mean number R of repetitions of hits on the same OM produced 
by direct and delayed photons in the detected muon-neutrino events, is used. The 90% c.l. upper 
limit for a E~ 2 energy spectrum is E 2 <3>go% — 5.3 x 10~ 8 GeV cirT 2 s _1 sr~' in the energy range 
20 TeV - 2.5 PeV. Other models predicting cosmic neutrino fluxes with a spectral shape different 
from E~ 2 are tested and some of them excluded at a 90% c.l.. 
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Model 


R* 




AEgo% 

(PeV) 




MPR[5] 


1.43 


3.0 


0.1-=- 10 


0.4 


P96py[30] 


1.43 


6.0 


0.2-r 10 


0.2 


S05[31] 


1.45 


1.3 


0.3-=- 5 


1.2 


SeSi[32] 


1.48 


2.7 


0.3-r 20 


0.6 


Mpp + py[33] 


1.48 


0.24 


0.8-=- 50 


6.8 



Table 3: Astrophysical flux models, the value of the R* which minimizes the MRF, the expected number of 
events N mod , the energy range AE 90 % in which the 90% of events are expected, and the ratio pt)Q%IN mol j. 
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